Breakdown of Conformal Invariance at Strongly Random Critical Points 
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We consider the breakdown of conformal and scale invariance in random systems with strongly ran- 
dom critical points. Extending previous results on one-dimensional systems, we provide an example 
of a three-dimensional system which has a strongly random critical point. The average correlation 
functions of this system demonstrate a breakdown of conformal invariance, while the typical correla- 
tion functions demonstrate a breakdown of scale invariance. The breakdown of conformal invariance 
is due to the vanishing of the correlation functions at the infinite disorder fixed point, causing the 
critical correlation functions to be controlled by a dangerously irrelevant operator describing the ap- 
proach to the fixed point. We relate the computation of average correlation functions to a problem 
of persistence in the RG flow. 



I. INTRODUCTION 

The principle of conformal invariance [jjj has proven 
immensely powerful in the study of continuous phase 
transitions. While it is believed to hold at critical points 
of systems with short ranged interactions in all dimen- 
sions ||], it is especially powerful in two dimensions, 
where it has enabled a substantial classification of crit- 
ical theories based on the representation theory of the 
Virasoro algebra and its extensions ||. 

Recent interest has focussed on applying this analy- 
sis to systems with quenched randomness in two dimen- 
sions; see for examples Refs. [||-^). The formal device 
that makes this possible is the construction of transla- 
tionally invariant field theories whose correlation func- 
tions are the disorder averaged correlators of the random 
problem, quite generally via the replica trick and in Gaus- 
sian problems, via the supersymmetry method. (There 
is at least one notable example where neither is needed 
0, and conformal invariance follows straightforwardly.) 

Our interest in this paper, is in asking whether the av- 
eraged field theories are necessarily conformally invariant 
at their critical points. We will find that this is not the 
case in two instances of what we term strongly random 
critical points, a category that may well have consider- 
able overlap with the "infinite disorder" critical points 
studied by means of real space renormalization group 
techniques [g]. Both are localization problems, the first a 
well studied problem in d = 1 where conformal invariance 
is even more powerful than in d = 2, and the second is 
an analogous construction on our part in d = 3. We also 
comment on another well studied member of this family 
of problems in d = 2 which does appear to have a region 
where conformal invariance holds. We should note before 
proceeding further, that while we were motivated by the 
absence of conformal invariance in the one dimensional 
example, what we find is a breakdown of scale invariance 
itself in that the operator product expansions is anoma- 
lous even though two point correlators are algebraic. 

Another question appears to be closely connected to 
the above issues: namely, how one might recover the dis- 



tributions for correlation functions that should charac- 
terize the universal content of random fixed points, from 
computations of averages and higher moments. In both 
of our examples, the typical correlations are very different 
from the average ones and are not even algebraic, which 
suggest both the difficulty of reconstructing them from 
the averaged field theory and why the latter may have to 
be anomalous. This raises the question of whether the 
breakdown of conformal invariance we report may have 
echoes at other, less strongly random critical points. We 
will speculate briefly on this at the end of the paper. 

One rationalization of our results will be that the av- 
eraged correlations vanish at the infinite disorder fixed 
point, causing the behavior of the critical correlation 
functions to be controlled by a dangerously irrelevant 
operator describing the approach to the fixed point. For 
the typical correlation functions, scale invariance breaks 
down, while the average correlation functions have power 
law behavior at the level of two point correlators, but 
break scale invariance for three and higher point correla- 
tions. 

Within this interpretation, the power law behavior of 
two-point correlation functions is not the result of a non- 
vanishing scaling dimension for the given operators at the 
critical point, as in this case one would find scaling for 
multi-point correlations as well. Instead, the power law 
behavior is a feature of the leading corrections to scaling. 

We begin by reviewing a one-dimensional system ex- 
hibiting the given breakdown in scale and conformal in- 
variance, and show the anomalous operator product ex- 
pansion. Next, we introduce a model system in three 
dimensions. We present evidence that the system is at a 
critical point, with relevant perturbations that introduce 
a localization length that diverges as the critical point 
is approached; however, we will be unable to describe 
fully the system off-criticality. We then compute typi- 
cal correlation functions in this system and show that 
they violate scale invariance. Using a Liouville field the- 
ory to compute average correlation functions, we show 
that while the two point functions are power law, indica- 
tive of a critical point and compatible with scale invari- 
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ance, multi-point correlation functions violate conformal 
invariance. 

In an Appendix, we note that there exists a SUSY 
field theory which reproduces the correlation functions 
of the Liouville field theory. Although the correlation 
functions are easier to compute in the Liouville theory, 
the advantage of the SUSY theory is that it provides 
a purely local field theory which exhibits a breakdown 
of conformal invariance; the Liouville field theory will 
involve one global integral which effectively introduces a 
long-range interaction. 



II. BREAKDOWN OF CONFORMAL 
INVARIANCE IN ONE DIMENSION 

The conformal group is the group of transformations 
which leaves angles unchanged. In three, or more, dimen- 
sions, the conformal group is generated by translations, 
rotations, dilatations, and inversions (x^ —> i§p)- By 
composing inversion-translation-inversion, one generates 
the so-called special conformal transformations. In two 
dimensions, the conformal group is supplemented by all 
analytic transformations of the complex plane, while in 
one dimension any diffeomorphism is a conformal trans- 
formation. 

Invariance under translation, dilatation, and rotation 
fixes the two point function to be a power law. Invariance 
under inversion fixes the form of the three-point correla- 
tion function. For example, consider a real scalar field, 
6, in d > 2 dimensions with action 



(1) 



Consider the connected two-point function of (f> 2 , 
(<f) 2 (0)(j) 2 (x)) cx |x| 4_2d , a power law as expected. The 
connected three-point function is given by 



(</> 2 (O)0 2 (^ 2 (y)) cx \x\ 2 - d \y\ 2 - d \x - y 



2-rf 



(2) 



in agreement with inversion symmetry. 

In one and two dimensions, conformal invariance im- 
poses even more stringent requirements on the correla- 
tion functions. In one dimension, conformal invariance 
requires that, at a critical point, all correlation functions 
are constant. An example of this behavior in statistical 
mechanics is the one-dimensional Yang-Lee edge || . 

Another example is the following quantum system. 
Consider a one-dimensional Dirac particle with Hamil- 
tonian 



H 



m(x) H 
m(x) — d x 



a, 



(3) 



If we take m(x) = m, then, looking at the zero energy 
Green's function, the point m = is a critical point. In 
accord with conformal invariance, the Green's function 



is a constant. For non-zero m, the system acquires a 
correlation length proportional to m _1 . 

Turning to a random version of this critical point, how- 
ever, we find a breakdown of conformal invariance. Let 
m(x) = m + 5m(x), with 5m{x) a quenched Gaussian 
random variable with vanishing mean and 



Sm(x)Sm(x') = g8(x — x' 



(4) 



where the overline denotes averaging over disorder. With 
this random mass, the Dirac equation is the continuum 
limit of a tight-binding Hamiltonian on a one-dimensional 
chain, with randomly chosen hopping elements and no 
on-site potential jl^] . The zero energy eigenfunctions of 
this problem provide an exception to exponential local- 
ization in one-dimension. 

For this problem, we will being by summarizing old 
results on the two-point functions. We will note that the 
two-point function, while invariant under translation, di- 
latation, and inversion, is not invariant under arbitrary 
diffeomorphisms, thus violating conformal invariance in 
one dimension. Then, we will show how the correlation 
functions may be obtained from the exact zero energy 
eigenfunctions of the problem, and use this to compute 
the three-point function. This function will not be invari- 
ant under inversion or special conformal transformations, 
thus breaking conformal invariance in a way that gener- 
alizes to higher dimensions. 

For m = 0, it has been shown [p^-|l5| that the average 
zero energy Green's function decays as a power law 



G(x,y) 



y\ 



-3/2 



(5) 



All positive moments of the average Green's function be- 
have similarly. For convenience, in this paper we choose 
to focus on a density of states correlation function. Let 
ipa(x) be an eigenfunction of Hamiltonian (||), with en- 
ergy E a . Introduce an infinitesimal parameter rj. Define 
the density of states p(x) to be equal to 



p(x) = Im(G(x,x)) — ImQj 



f 



E„ 



IT] 



:\Mx)\ 2 )- (6) 



The two point correlation function of the density of 
states scales as Eq. (^). On the other hand, if we look at 
the full probability distribution of the density of states 
at two points, we find that if p(x) is of order unity, then 
the typical p(y) is of order 



-\fg\x- 



(7) 



Eq. (Q) is not invariant under general diffeomorphisms, 
while Eq. ([?]) is not scale invariant. However, we empha- 
size that the point m = is a critical point for the random 
system, as for m 7^ 0, the average Green's function ac- 
quires a diverging correlation length proportional to m -2 , 
while the typical correlation function is proportional to 
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to . Similarly, away from zero energy one finds a finite 
localization length, diverging as ln 2 (l/£:) or ln(l/£T) for 
average and typical correlation functions, respectively. 

To derive the above results, consider the exact zero 
energy eigcnfunctions of the problem: 







(8) 



where f3{x) = J m(y) dy. Representing the zero energy 
Green's function as a resolvent 



G(x,y) = (x\ 



1 



H + irj 



(9) 



we note that, for a finite system in the limit r] — > 0, 
the Green's function and density of states are controlled 
by the zero energy eigenfunction; the imaginary part of 
the Green's function is equal to -ip(x)i/j(y) where tp is 
the zero energy eigenfunction. The zero energy Green's 
function is a matrix; the different components of this ma- 
trix provide information on the two different zero energy 
eigenf unctions. In the Appendix, a supersymmetric field 
theory is constructed such that the zero energy Green's 
function is a correlation function of the form 



(ip(x)ip(y)). 



(10) 



If the infinite volume limit is taken before the i] — > 
limit, the resolvent and density of states depend on other 
eigenfunctions beyond the exact zero energy eigenfunc- 
tion. We now show, however, that the scaling properties 
of correlation functions for both these quantities are the 
same for either order of limits. Consider eigenfunction 
ip a at energy E a . The eigenfunction is a zero energy 
eigenfunction of Hamiltonian 



H = 



-E m(x) + d x 
m(x) — d x —E 



By performing an axial transformation, 

• e -m o 



fpa 







o/3(s) 



(11) 



(12) 



we find that ip a is a zero energy eigenfunction of Hamil- 
tonian 



H 



_ Ee 2f3(x) Q x 

-d x -Ee~ 2 ^ 



(13) 



This Hamiltonian has diagonal disorder, and one expects 
to find all the eigenfunctions localized; however, for small 
E, the diagonal terms can be ignored for lengths such 
that Ee 2/3 << 1. Since (3 has rms fluctuations of or- 
der the diagonal terms in Hamiltonian ( p"3| ) can be 
ignored for L « \ii 2 (l/E), and the exponential factor 
from the axial transformation (^) controls the magni- 
tude of ip a (we note that the exponential behavior of 



P ensures that the dominant contribution to the mag- 
nitude of ip a arises from the axial transformation, rather 
than any fluctuations in magnitude of ip a ). Thus, we find 
that over short length scales all the eigenfunctions track 
the zero energy eigenfunction in relative magnitude; over 
longer length scales, the non-zero energy eigenfunctions 
can become localized. 

Now consider correlation functions of density of states 
for finite to Since the density of states includes eigen- 
functions with energy E of order to we find that rj sets a 
length scale in the system (this length is identical to the 
average localization length discussed above). However, 
on lengths shorter than this scale, the magnitude of of the 
eigenfunctions contributing to the density of states can 
be obtained from that of the zero energy eigenfunction, 
up to constant factors. Consider a correlation function of 
the density of states, [j/>(ii). We argue that these cor- 

i 

relation functions, in the thermodynamic limit with fixed 
to are equivalent, up to constant factors, to correlation 
functions taken in a finite system as rj — > 0, provided that 
in the first case \xi — xj\ << ln 2 (i). A similar argument 
may be made for average Green's functions. 

Given this equivalence, for the remainder of the paper 
we consider only the behavior of the zero-energy eigen- 
function. The log of the exact zero-energy eigenfunction, 
(3{x), executes a random walk, such that the mean-square 
fluctuations in the log are proportional to the length scale 
L. If we focus on the first of the two eigenfunctions, we 
see that it is strongly peaked at a given point, the max- 
imum of (3{x) (the other eigenfunction is peaked at the 
maximum of — (3(x)). 

After normalizing the eigenfunction, it is exponentially 
small away from the maximum. Typically the eigen- 
function decays as a stretched exponential, giving Eq. 
(0). It is possible, though, for the eigenfunction to have 
a secondary maximum. This possibility dominates the 
average two-point function. Eq. (^|) can be obtained 
as follows: we wish to compute the probability that 
(3{x) has global maxima at 0, x. Since the eigenfunction 
normalization can be absorbed by a shift in /3, we fix 
/3(0) = (3{x) = 0. Everywhere else, (3{y) < 0. Thus, the 
probability to have maxima at 0, x scales as the return 
probability of a random walk required to persist below 
0. This scales as |x|~ 3 / 2 , giving Eq. (||). Further, the 
random walk technique shows that all positive moments 
have the same scaling, as any positive moment of the zero 
energy eigenfunction is sharply peaked at the maxima of 
p. 

Tuning away from the critical point by taking m ^ 0, 
we find that the zero energy eigenfunction is normaliz- 
able only for a finite system, where it describes eigen- 
functions decaying into the system. The decay of the 
eigenfunctions is governed by the diverging correlation 
length discussed above. 

We now turn to the three-point correlation function 
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of the density. This function can be represented in the 
supersymmetric field theory as 



(14) 



The average three-point correlation function can be com- 
puted by an extension of the random walk argument. As- 
sume that < x < y. Then, we consider a random walk 
that starts at 0, returns to zero at x, and returns again 
to at y. This gives the average three-point correlation 
function 



p(0)p(x)p(y)<x\x\- 3 / 2 \x-y\- 3 / 2 . 



(15) 



As stated above, Eq. Pq ) is not invariant under inversion 
or special conformal transformations. 

The three-point function enables us to find the oper- 
ator product expansion for the theory. Given a set of 
operators, Oa, with dimensions, cLa, one expects that 



O A (0)O B (x) -> fABcx- [dA+dB - dc) Oc 



(16) 



Since all positive moments of the density of states scale 
equivalently, we need consider only two operators. One 
operator is the identity operator, /, while the other oper- 
ator O is the density of states operator, with dimension 
d. Eq. (||) suggests d = —3/4. Then, consistency with 
the operator product expansion would give 

P(0)p(x)p(y) ex |*r 3 / 4 |* - y|- 3/4 |yr 3/4 , (17) 

which, like Eq. (^|), obeys inversion symmetry. How- 
ever, this is not what we have found. Thus, we ratio- 
nalize our results as follows: the scaling dimension of 
the operator O vanishes, as indeed it must in one di- 
mension. However, the rms fluctuations in (5 scale as 
and the inverse rms fluctuations in (3 tend to zero 
as L — ► oo; equivalently, the statistical properties of the 
zero energy eigenfunction are unchanged under scaling 
x — > x/L,g — ► Lg. If, however, we take g divergent, we 
find that the zero-energy eigenfunction is localized at a 
point and all correlation functions vanish. Thus, g~ l is a 
dangerously irrelevant variable. In the scaling limit, de- 
fined as the limit in which the separation between points 
in correlation functions tends to infinity, all correlation 
functions are controlled by the vanishing of this variable 
as the fixed point is approached. Similarly, /abc van- 
ishes as the fluctuations in j3 diverge. However, what we 
have found above is that / vanishes as a power law as 
the fixed point is approached. Thus, the power law (Q) is 
not the result of d + 1 0, but the result of / — > as the 
fixed point is approached. 

We note that the behavior we have found is similar to 
that found using a real-space RG approach []l3| on the 
discrete tight-binding version of this model. In that ap- 
proach, one renormalizes the distribution of hopping ele- 
ments, reducing a cutoff on the strength of the hopping. 



In the scaling limit, almost all the hopping elements be- 
come vanishingly small compared to the cutoff; the frac- 
tion of hopping elements (vanishing as a power law in 
the scaling limit) that are of order the cutoff control the 
correlation functions. 



III. A MODEL SYSTEM IN THREE 
DIMENSIONS 

To demonstrate the generality of this breakdown of 
conformal invariance, we will look for higher dimensional 
examples. Using a real-space renormalization group, 
some examples have been found of infinite disorder quan- 
tum critical points two or more dimensions flq| . How- 
ever, these are interacting quantum systems and can only 
be described by 2+1 or 3+1 dimensional theories with a 
privileged time coordinate, so they would not be expected 
to have conformal invariance, although the question of 
scale invariance in these systems remains an interesting 
problem for future work. There is also some evidence 
|]l7f that a two-dimensional random hopping model could 
have an infinite disorder critical point. However, there is 
as yet no evidence for breakdown of conformal invariance 
in this model. Therefore, we will construct a model sys- 
tem in three-dimensions that provides a strongly random 
critical point without conformal invariance. 

After presenting the Hamiltonian and eigenfunctions, 
we will present evidence that the system is indeed at a 
critical point between two different localized phases. Fur- 
ther evidence for the criticality of our model will result 
from behavior of typical and average correlation func- 
tions that we find below. 

Consider the following Hamiltonian 



H 



where 



(0„ 



A fl {x)=d fi /3(x). 



A,) 2 



(18) 



(19) 



The Hamiltonian has two exact zero energy eigenfunc- 
tions. They are 











(20) 



We pick (3{x) with probability distribution e s ^ 
where 



S\P] = Yg / (,Kri " 2 " 2 



(21) 



We pick this distribution of disorder as the best balance 
between having too little and too much disorder, while 
retaining a local distribution of disorder. Let us show 
that for less disorder, the critical point would not be 
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strongly random, while for more disorder we might lose 
the critical point completely. 

With this distribution for /J, the mean-square fluctua- 
tions in /3 are proportional to L, the length scale of the 
system, just as in the one-dimensional system the fluctua- 
tions in J m(y) dy were proportional to the length of the 
system. This increase in fluctuations with length scale 
enables us to reach a strongly random fixed point. 

Suppose instead that we had chosen a disorder distri- 
bution with smaller fluctuations in (3. If we had cho- 
sen a two-dimensional system with S[(3] = J cPx (d^fi) 2 , 
the fluctuations in (3 would only grow logarithmically 
with L, the exponential of [3 would behave as a power 
law, and we find a conformally invariant fixed point [fl8f . 
If we had considered a three-dimensional system with 
S[j3] = J d 3 x (d^P) 2 , we would have found a fixed point 
with vanishing disorder. 

On the other hand, for the given distribution of dis- 
order, we can make an argument that the system is at 
a critical point. For the given distribution of disorder, 
by letting = d t _ l /3(x) + J M , for J M ^ 0, we can tunc 
away from criticality. As in the one-dimensional exam- 
ple, we find that the eigenfunction is normalizable only 
for a finite system and describes eigenfunctions decaying 
into the system. As now the zero energy eigenfunctions 
are 









-0{x)- 



JuX" 



(22) 



the typical eigenfunction decays with a length propor- 
tional to |J| _1 - This demonstrates that there is a di- 
verging length scale associated with the given problem, 
implying that we are at a critical point. If we consider 
a system in a square box, with = (J, 0, 0), as we tune 
j from positive to negative there is a phase transition 
between eigenfunctions localized on opposite faces of the 
box. It should also be possible to tune away from critical- 
ity by considering Green's functions at non-zero energy. 
This is a more difficult task, for future work. 

However, if we had considered a two-dimensional sys- 
tem with S[0\ — J d 2 x {d 2 f3) 2 , the mean-square fluctu- 
ations in f3 would grow as L 2 . We would find that the 
typical correlation functions for eigenfunction ( po| ) would 
decay exponentially. In the case of the system in a box, 
even for vanishing J there would be eigenfunctions ex- 
ponentially localized on the faces of the box, and the 
phase transition would occur at a non- vanishing value of 
J. Therefore, for this two-dimensional distribution we 
would not have such strong evidence that the system is 
critical. 

Finally, let us emphasize that the distribution of dis- 
order we have chosen in Eq. (|2l]) is local. As a result, 
the SUSY field theory in the Appendix is also local, and 
so the breakdown of conformal invariance is not a result 
of long-range interactions. 



IV. TYPICAL CORRELATION FUNCTIONS 
AND ONE-POINT FUNCTIONS 

Let us focus on the first zero energy eigenfunction of 
Eq. rt2Q) and write ip(x) = e 13 ^. We normalize tp such 



that J d 3 xip 2 (x) = 1. In exact analogy with the one- 
dimensional system, we find that if tp(x) is of order unity, 
then 



(23) 



We can also consider a typical one-point function, the 
magnitude |i/'(a;)| 2 . Of course, on average, this function 
is equal to -p-. Typically, however, it is of order 



(24) 



Clearly, Eqs. (|23j,|24j) violate scale invariance. We will 
next turn to the average correlation functions. They will 
not violate scale invariance, but they will violate con- 
formal invariance. The advantage of considering average 
correlation functions is that, as shown in the Appendix, 
they can be obtained as correlation functions in a purely 
local field theory. 



V. MULTIFRACTALS AND LIOUVILLE FIELD 
THEORY 

Let us look at average correlation functions of the form 

(25) 



W q (x,L) = \iP(x)iP(0)\*, 



where the overline denotes ensemble averaging. From 
the scaling behavior of W q , we can compute the scaling 
of p{x)i/ 2 p(0)i/ 2 . 

The inverse participation ratio, W q (L) — W q (0,L), is 
expected to vary as a power law 



W q (L) cx 



f 



L 3+r(q) ' 

while the two-point function is expected to vary as 

W„(x,L) OC -z , , -r-r. 



(26) 



(27) 



For localized eigenfunctions, r{q) — 0, while for a plane 
wave eigenfunction r(q) = 3(q — 1). In some localization 
problems, one finds multifractal scaling for the r(q) at 
criticality; see ref. Jl!| for a review. 

The typical correlation functions of the previous sec- 
tion can be considered as exponentials of ensemble av- 
erages of logarithms of correlation functions. They can 
be obtained from q — > limits of the average correlation 
functions. 

We will now introduce a Liouville field theory [pO 14 
to compute the functions W q . The correct normalization 
of eigenfunction ip can be obtained by 
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^{x)=N- l ' 2 e^ x \ (28) 



N = / d 3 xe 2l3( - x \ (29) 



where 



The correlation functions W q can be written as 

W q {x,L) = ±r [[dflN-^WeiK'h-sW, (30) 
Zo J 

where Zq = f[dp]e~ s ^ . The normalization can be ex- 
ponentiated by 



N~ q 



r(«) 



duu t ~ 1 e- uN . 



(31) 



As a result, 



Wq (x,L) = ( dojm^e^e^e- 3 ^-^ 

A)l (?) 



2/3 







(32) 



This provides a three-dimensional Liouville field theory, 
with the strange feature that the action involves fourth 
derivative terms instead of second derivative terms. 



VI. AVERAGE CORRELATION FUNCTIONS 

We will consider the average correlation functions. We 
will demonstrate that the exponents r(g) and <j{q) are 
independent of q. We will be unable to calculate a(q) 
exactly, but we will argue that the power law behavior of 
Eq. ( p7| ) is obeyed for a finite cx(q), leading to conformal 
invariance for the two point function. The behavior of the 
average correlation function is again related to problems 
of persistence in noncquilibrium systems. Then, as in one 
dimension, we will show using OPEs of the theory that 
the three point function breaks conformal invariance. 

If we look at the scaling behavior of Eq. (|32"1), we find 
that as the length scale L increases, (3 increases propor- 
tional to L 1 / 2 , and so e q @ becomes sharper. Effectively, q 
increases proportional to L 1 / 2 , and so for any q > 0, the 
scaling behavior of the correlation functions is the same, 
so that r(q) and cr(q) are independent of q. The q — > 
and L — ► oo limits do not commute. 

In the large q limit, and therefore in the scaling limit, 
the two-point function is dominated by field configu- 
rations in which (3 reaches its maximum at and x. 
The one-point function is dominated by configurations 
in which (3 reaches its maximum at 0; the probability of 
this happening is proportional to 1/L 3 and so r(q) = 0, 
characteristic of a localized eigenfunction. 



To find cr(q), we need to look at the scaling of the 
following functional integral 



-sw-flP) 



(33) 



Here, Pl [j3] = 5((3(0)) and p 2 [[3] = S{f3{x)), while f{J3) = 
for (3 < and f((3) = oo for (3 > 0, so that the potential 
/ has a hard wall [0 at (3 = 0. This reproduces the 
desired constraint that (3 have a maximum at both and 
x. Let us define an operator O(x) that inserts a factor of 
6(f3(x)) in Eq. @. Then W q (x,L) = (O(O)O(x)). 

We will proceed with a renormalization group tech- 
nique as follows: initially the functionals pi,p2 define 
(5-function distributions for (3 at 0, x. We will replace 
these with probability distribution functionals for [3, and 
for the normal derivative of (3, on the surface of a sphere 
at radius a around points 0, x. Initially a will be small, 
but we will integrate out the field j3 in a shell between 
a and a + da, and compute the change in the probabil- 
ity distribution as a is increased. In one-dimension, this 
amounts to computing the probability distribution of a 
random walk, persisting below zero, at points 0±a, x±a. 

As the RG above is carried out with spheres centered 
on cc,0, at a = \x/2\ the two spheres merge, and one 
must match boundary conditions. On the surface of the 
merged spheres, the field (3 is of order To compute 

the scaling of the two point function, one must compute 
the probability that (3 persists below zero on the surface 
of both spheres up to scale |x/2|, with matched bound- 
ary conditions (we will discuss this computation below). 
After the spheres merge, one has a single surface prop- 
agating outwards as a is increased, but the persistence 
probability on this surface is independent of x. 

Given the two-point function, consider the three-point 
correlation function (O(0)O(x)O(y)), with << \y\. 
At the scale \y\, the field (3 is of order >> \/\x\- 

Therefore, for the purpose of computing the scaling of 
the three-point function with respect to y, we can simply 
assume that (3 is pinned at the value zero near the point 
0. Thus, the two operators O(0), 0(x) fuse into a single 
operator O. 



lim O(x)O(0) = x- 

x—>-Q 



r O(0). 



(34) 



The exponent as a in the above equation is exactly cr{q) 
for q > 0, as the two-point function can be computed by 
fusing the two operators O into a single operator, yielding 
a one-point function, for which 



Therefore, 



(O(0)) = 1/L 3 



(O(O)O(x)) cx x- 



(35) 



(36) 



We now discuss the procedure for computing a. First 
we consider the problem of persistence in (3 for a single 



G 



sphere, and then show how to obtain a from the per- 
sistence problem by matching boundaries. However, the 
problem in persistence is sufficiently difficult that we will 
not attempt to obtain a numerical value for a. The re- 
lated persistence problem will at least enable us to argue 
that a is finite so that correlation functions are power 
law. 

Decomposing [3 into spherical harmonics, 



p(a,e,4>)=Y J fim{a)Y lm ( 



(37) 



l,m 



the probability distribution functionals can be written as 



p[fim(a),d a fi m (a)}. 



(38) 



First, let us consider the probability that, for a single 
sphere, j3 does not cross the hard wall up to scale a. 
Defining 

9im(T) = e- T ' 2 f lm {e T ), (39) 
we can write the action S[(3] as 

Y, J ^(3/4 + Z(Z + l) + 2d T + c>|) ff/m (T)) 2 dT. (40) 

Ira 

With this scaling of field, and logarithmic transform of 
the RG scale, the problem becomes that of persistence in 
a Gaussian stationary process (2^] . With action given by 
Eq. (pi), the correlation function 



(g lm (T)g lm (T')) 



(41) 



is exponentially decaying in |T — T'\. Therefore, for a 
finite number of fields g 7 the persistence probability is 
exponentially decaying in T, and decays as a~ s for some 
power 8. We have an infinite number of fields gi rrll but 
for large I they have very small fluctuations, so it seems 
likely that 8 tends to a finite limit as we include fields 
with greater and greater I. 

In principle it is possible to compute 8 for a finite num- 
ber of g. However, since the process is non-Markovian 
due to the fourth derivatives it is quite difficult; even 
for the simpler problem of one field, with an action con- 
taining only fourth derivative and no other terms, this 
problem [^3| was solved only relatively recently jgjf ]. 

For a finite system, without periodic boundary condi- 
tions, the persistence exponent determines the probabil- 
ity that the maximum of the eigenfunction will be located 
at a given interior point, instead of at the boundary. 

Given the persistence exponent, we can compute a. At 
scale a = \x/2\, the spheres start to merge. Beyond this 
scale, we have only one surface. Immediately before the 
spheres touch, the probability that (3 is a maximum at 
in the sphere centered around is 

\x\- e . (42) 



The probability that f3 is a maximum at x in the sphere 
centered at x is the same. The product is |cc| -2e . 

However, at the point where the two spheres meet, the 
values of (3 and d(3 must be equal on both spheres. This is 
a matching of boundary conditions. To match boundary 
conditions, the scaled variables g must be equal on both 
spheres. Since g is stationary, the probability that the 
scaled variables arc equal is independent of the RG scale 
T. However, for /3 to be equal on both spheres, there is an 
additional Jacobian which gives rise to a further power 
law correction. In this case, the mean square fluctuations 
in j3 grow as \/L, while the fluctuations in d(3 do not grow 
in L. Therefore, the probability that (3 is a maximum at 
both and x, and that the boundary conditions match 
(equivalently, that the maxima at and x coincide) is 

This is to be compared to Eq. ( ff2| ) for only one sphere, 
so that the relative probability to find two maxima is 



a = 28 + 1/2 



1/2. 



(43) 



We see this explicitly in the problem of the random Dirac 
equation (^) in one-dimension, where the persistence ex- 
ponent 8 = 1 (typically in the literature, this is quoted as 
8 = 1/2 as one typically computes only the persistence 
probability for a noncquilibrium process going forward in 
time; our spheres move outwards in both directions and 
so the exponent is doubled). The average correlation 
function decays as a;~ 3 / 2 = x~ e ~ 1 / 2 . 

Even without a precise value for cr, we can still demon- 
strate the breakdown of conformal invariance. The two- 
point function calculated above exhibits scale and con- 
formal invariance in 3 dimensions. However, as in one 
dimension, the operator product expansion is not con- 
sistent with scaling. From the two-point function, one 
would expect that O would have scaling dimension a/2. 
However, Eq. ( |34| ) gives a scaling dimension of a for O. 
The trouble is that Eqs. (|i 36) are not simultaneously 
consistent with scaling. 

To make the contradiction more concrete, consider the 
three point function W q (x,y.L) = \tp(0)ip(x)ip(y)\i = 
(O(0)O(x)O(y)) in the limit when \x\ « \y\. Although 
we were unable to compute cr, we are able to use the 
operator product expansion to demonstrate a breakdown 
of conformal invariance. 

In the limit \x\ << \y\, the operator product expansion 
gives 



W q (x iy ,L)KL- 3 \x\-°\y\- 



(44) 



This is similar to Eq. (|15|) found in one dimension. Eq. 
( p4| ) violates inversion symmetry with respect to a point 
\z\ where \z\ « \x\. 

As in one-dimension, we interpret this to mean that O 
has vanishing scaling dimension, but that all the corre- 
lation functions of O vanish at the infinite disorder fixed 
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point. As the theory approaches the fixed point, the av- 
erage correlation functions vanish as a power law, so that 
if one only considers the two-point functions one obtains 
the wrong scaling dimension for O. The inverse of the 
mean-square fluctuations in [3 is a dangerously irrelevant 
operator that describes the approach to the fixed point; 
when this operator vanishes, so do the correlators. 

VII. CONCLUSION 



only contribution to an n-point correlation function of 
spins will arise when all n points are in a locally fer- 
romagnetic region. In this case, we conjecture that the 
computation of correlation functions will again become a 
problem of persistence in the RG flow, that the system 
will stay within the ferromagnetic phase. Therefore, we 
conjecture that a breakdown of conformal invariance will 
be seen for sufficiently large moments in other systems 
as well. 



We have considered the breakdown of conformal in- 
variance in certain strongly random fixed points. The 
breakdown of conformal invariance is related to a power 
law vanishing of correlation functions approaching the 
fixed point. 

Although we were unable to compute the exponent 
a precisely, we were able to obtain some understand- 
ing of the operator product expansion, and to relate a 
to a persistence exponent in a noncquilibrium process. 
The connection to persistence is not surprising; the RG 
flow describing the approach to the fixed point is non- 
equilibrium, as the flow approaches the fixed point but 
never reaches it. In contrast to ordinary critical points, it 
is the approach to the fixed point, rather than the fixed 
point itself, that controls the exponents. 

Within other strong disorder RGs, such as for the 
transverse field Ising model |I^,|l6|], average correlation 
functions are controlled by the probability that a given 
pair of sites will remain within the ordered phase as the 
RG is run up to the scale of the separation between the 
sites. This is also a question of persistence. 

Interestingly, for other random fixed points, such as 
the problem of two-dimensional Dirac fermions in ran- 
dom magnetic field jl8| or the problem of Anderson lo- 
calization it has been found that averages of low 
moments of correlation functions exhibit conformal in- 
variance, but for sufficiently high moments conformal in- 
variance breaks down. Shapiro showed for the Anderson 
localization problem [[ll| that the breakdown of confor- 
mal invariance could be described by a universal distri- 
bution of conductance with a broad tail, such that higher 
moments were divergent when evaluated with the fixed 
point distribution of conductance. As a result, the higher 
moments are controlled by the approach to the universal 
distribution. 

We would like to close with a conjecture. For the 
strongly random systems we consider, all positive mo- 
ments of the correlation functions have the same scal- 
ing behavior. For systems with weaker randomness, one 
observes different scaling dimensions for the various mo- 
ments, q. However, there is an expectation that the scal- 
ing dimension will saturate, and become independent of 
q for sufficiently large q. At this point, one may cal- 
culate correlation functions with q = oo. Consider for 
dcfiniteness a disordered Ising model. In this case, the 



VIII. APPENDIX: A SUSY FIELD THEORY 

In this Appendix, we will show that there exists a su- 
persymmetric field theory 1 21 1 which reproduces the aver- 
age correlation functions computed above in the Liouville 
field theory. Variants of the supersymmetric field theory 
can be used to describe the one-dimensional and three- 
dimensional problems considered. The construction of 
the supersymmetric field theory will show that there ex- 
ists a purely local field theory with a finite number of 
fields which violates conformal invariance at the critical 
point. In the Liouville theory, the integral over to in 
Eq. ( |32] ) introduces a non-local interaction, which might 
otherwise appear to be responsible for the breakdown of 
conformal invariance. 

We are interested in the resolvent 



G{x,y) = (x\ 



1 



H 



irj 



:\y)- 



(45) 



If the 77 — ^ limit is taken before the infinite volume 
limit, the resolvent is determined by the exact zero- 
energy eigenfunction. As discussed above, all the local- 
ization properties near zero energy are controlled by the 
zero energy eigenfunction. 

The computation of averages of G follows standard 
techniques f2l]| . Introducing a superfield with bosonic 
component <p and fcrmionic component rp, the resolvent 



is 



where 



e -S[/3] e -*(iH-n)« 



Z 



(46) 



(47) 



By introducing a set of m superfields <&i, we can com- 
pute m-th moments of the correlation functions. Let us 
emphasize that we can take m to be fixed, but large, 
and then we have a single field theory from which high 
moments and multi-point correlation functions can be 
computed. 

In order to compute the typical correlation functions of 
section IV, we would need to introduce an infinite set of 
superfields to compute all moments, and then obtain the 
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distribution of correlation functions from the moments. 
Alternately, we could compute them from the Liouville 
field theory as a limit of the q-th moment of the correla- 
tion functions as q — > 0. In any event, we do not know 
how to obtain them from a purely local field theory. 

In the SUSY field theory, (3 does not acquire any dy- 
namics as $ is integrated out. Physically, this is the 
statement that (3 is quenched. If (3 were not quenched, 
one would expect (3 to acquire a second derivative term 
under the RG, and then the theory would flow to a fixed 
point with finite fluctuations in (3. The fact that disor- 
der is quenched is important for the existence of a strong 
disorder critical point. 
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